时 变 能 量 网 络 理论 


陈 卑 勇 ， 葛 海 麟 ， 名 明 
《华南 理工 大 学 电力 学 院 广州 510641) 


摘要 


不 同类 型 的 能 源 系 统 通 过 能 量 转 换 设备 (感应 电动 机 、 离 心 泵 等 ) 相互 耦合 ， 研 究 多 类 型 能 源 系 统 的 
动态 特性 及 其 仿真 方法 对 多 能 互补 系统 优化 设计 及 性 能 分 析 具 有 重要 的 实际 意义 。 为 对 时 变 能 量 网 络 进行 
建 模 和 分 析 ， 本 文 从 能 量 本 质 的 角度 出 发 ， 通 过 深入 探讨 能 量 传递 及 转换 机 理 ， 分 别 建立 时 变 传递 线 〈 管 ) 
路 和 能 量 转换 设备 的 集中 参数 模型 。 在 时 变 能 量 网 络 模型 的 基础 上 ， 提 出 通过 构建 时 变 能 量 网 络 方程 ( 包 
括 状 态 方程 和 输出 方程 》 对 多 能 互补 系统 的 动态 特性 进行 建 模 仿真 的 分 析 方 法 ， 最 后 通过 具体 算 例 对 本 文 
所 提 分 析 方 法 的 有 效 性 及 实用 性 进行 验证 。 本 文 的 研究 内 容 为 时 变 能 量 网 络 的 建 模 、 分 析 、 优 化 及 规划 芮 
定 了 基础 。 
关键 词 : 广 延 量 ， 能 量 ;， 烟 ; 时 变 ， 能 量 网 络 方程 
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Time-varying Energy Network Theory 
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ABSTRACT: Different types of energy systems are coupled by energy conversion devices (e.g., induction motors, 
centrifugal pumps, etc.). It is of great practical significance to study the dynamic characteristics and simulation 
method of multi-type energy system for the optimization design and performance analysis of multi-energy 
complementary system. In order to model and analyze the time-varying energy network, from the point of view of 
the essence of energy, lumped parameter models of time-varying transfer line (pipe) and energy conversion 
equipment are established through in-depth study of the mechanism of energy transfer and conversion. On the 
basis of the time-varying energy network model, a method of modeling and simulating the dynamic characteristics 
of the multi-energy complementary system by constructing time-varying energy network equations (including 
state equations and output equations) is proposed, then the validity and practicability of which are verified by a 
practical example. The research of this paper lays the foundation for the modeling, analysis, optimization and 
planning of time-varying energy network. 
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1 引言 


实现 不 同形 式 能 源 的 综合 利用 ， 在 全 球 能 源 及 环境 危机 的 背景 下 ， 已 成 为 了 必然 的 趋 
势 " ”1 。 工 程 中 能 源 的 生产 、 传 输 及 利用 ， 大 多 采用 网 络 化 的 方式 进行 ， 如 电网 、 热 网 、 燃 和气 
管 网 等 。 英 国 和 爱尔兰 成 立 的 能 源 网 络 协会 ， 对 电力 和 燃气 输 配 管 网 的 协调 管理 展开 研究 D。 
瑞士 联邦 理工 学 院 在 “未 来 能 源 网 络 远 景 ” 项 目 中 ， 首 次 提出 了 Energy Hub 和 Energy 
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Interconnector 的 概念 1。 文献 [8] 指 出 由 可 再 生 能 源 、 电 力 系 统 和 信息 通信 技术 组 成 的 能 源 互 
联网 将 成 为 第 三 次 工业 革命 的 核心 。 能 源 互 联网 不 仅 实 现 电 能 生产 者 和 消费 者 之 间 的 共享 ， 
也 实现 不 同形 式 能 源 之 间 的 互联 和 共享 。 

不 同类 型 的 能 源 从 物质 的 角度 上 是 不 能 转化 的 “如 电 不 能 转化 为 落 汽 等 )， 然 而 从 能 量 的 
角度 却 是 可 以 相互 转换 的 (如 电能 转换 为 热能 等 ); 能 量 是 各 种 不 同形 式 能 源 与 物理 过 程 实现 
相互 转换 的 桥梁 1。 要 实现 不 同形 式 能 源 物 理 层面 上 的 互联 ， 必须 从 能 量 的 角度 出 发 , 探讨 
能 量 的 普遍 化 传递 规律 ,文献 [11] 在 深入 探讨 能 量 本 质 的 基础 上 , 建立 了 能 量 网 络 的 基本 理论 。 
然而 ， 工 程 上 能 源 系统 的 状态 往往 与 运行 工 况 有 关 ， 系 统 结构 或 运行 方式 的 改变 ， 都 会 使 系 
统 参数 动态 变化 ， 此 时 系统 参数 往往 与 时 间 相关 ， 因 此 非常 有 必要 进一步 建立 时 变 能 量 网 络 
的 相关 理论 。 

与 时 不 变相 比 ， 时 变 能 量 网 络 的 传递 过 程 将 出 现 一 些 新 的 规律 和 特点 。 本 文 从 能 量 本 质 
的 角度 出 发 ， 对 时 变 传递 线 〈 管 ) 路 的 普遍 化 传递 规律 进行 深入 探讨 ， 引 出 时 变 能 量 网 络 基 
本 参数 的 定义 及 物理 意义 。 结 合 不 同类 型 能 量 的 具体 传递 特性 ， 通 过 严谨 的 数学 推导 得 到 时 
变 传递 线 〈 管 ) 路 的 等 效 传递 方程 ， 建 立时 变 传 递 线 〈 管 ) 路 的 集中 参数 模型 。 不 同类 型 的 
能 源 系 统 通过 能 量 转换 设备 〈 感 应 电动 机 、 离 心 泵 等 ) 相互 耦合 ， 本 文 进一步 建立 能 量 转换 
设备 的 集中 参数 模型 。 为 考察 时 变 能 量 网 络 的 动态 过 程 ， 本文 在 时 变 能 量 网 络 模型 的 基础 上 ， 
提出 通过 构建 时 变 能 量 网 络 方程 《包括 状态 方程 和 输出 方程 )》 对 多 能 互补 系统 的 动态 特性 进 
行 建 模仿 真 的 分 析 方法 。 本 文 最 后 以 一 个 数值 算 例 证 明 所 建立 的 理论 及 方法 的 正确 性 及 实用 
性 ， 并 阔 述 了 其 基本 应 用 。 


2 时 变 能 量 网 络 的 基本 参数 


描述 热力 学 体系 的 状态 参量 可 以 分 为 强度 量 和 广 延 量 两 类 : 强度 量 是 指 与 物质 的 量 无 关 ， 
不 具有 加 和 性 的 量 ( 如 热力 学 温度 、 电 势 、 流 体 压强 ); 广 延 量 是 指 与 物质 的 量 有 关 ， 具 有 加 
和 性 的 量 〈 如 热力 学 焙 、 电 荷 量 、 流 体 体 积 )。 广 延 量 可 分 为 基本 广 延 量 、 能 量 及 烟 〈 如 电荷 
量 、 电 能 及 电 烟 ， 又 如 热力 学 粮 、 热 能 及 热 烟 )， 基 本 广 延 量 是 另外 两 种 广 延 量 的 传递 载体 。 
自然 界 任何 形式 的 能 量 和 烟 都 可 以 表示 成 一 个 基本 强度 量 和 一 个 基本 广 延 量 的 函数 ， 而 且 这 
两 个 基本 参量 的 乘积 正好 具有 能 量 的 量 纲 ( 如 电能 可 以 表示 成 电势 与 电荷 量 的 乘积 ， 热 能 可 
以 表示 成 温度 与 米 的 乘积 ) "|。 

为 进一步 引出 时 变 能 量 网 络 基本 参数 ， 下 文 将 从 基本 广 延 量 和 能 量 的 空间 平衡 方程 出 发 ， 
通过 严谨 的 数学 推导 ， 对 时 变 传递 线 ( 管 ) 路 的 传递 规律 进行 分 析 及 讨论 。 
2.1 时 变 传递 线 〈 管 ) 路 的 传递 规律 

能 量 网 络 由 传递 线 ( 管 ) 路 组 成 ， 因 此 须 准确 地 把 握 能 量 在 时 变 线 〈 管 ) 路 中 的 传递 规 
律 。 工 程 中 的 线 《〈 管 ) 路 通常 是 圆柱 形状 的 ， 如 电线 、 供 热管 道 及 天 然 气管 道 等 ， 本 文 将 以 
图 1 所 示 的 圆柱 形 传递 线 〈 管 ) 路 为 对 象 ， 对 时 变 传递 线 〈 管 ) 路 的 传递 规律 进行 深入 探讨 。 


A 
ea 
X X， 
图 1 圆柱 形 传 递 线 〈 管 ) 路 
Fig.1 Cylindrical transfer line (pipe) 


基本 广 延 量 和 能 量 的 空间 平衡 方程 分 别 为 021 
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dx 


Ty 1 

P J +g, (1) 
Xd. 

p= VT) + go O) 


其 中 jp 为 介质 的 密度 ，x 为 单位 质量 介质 中 的 基本 广 延 量 ，X 为 强度 量 ，J, 为 基本 广 延 


量 的 流 密度 矢量 ， g, 为 单位 体积 介质 中 基本 广 延 量 的 源 强度 ， g。,=J, .VX + Xg, 为 单位 体积 介 
质 中 能 量 的 源 强 度 。 对 上 两 式 分 别 进行 体积 积分 ， 则 有 


CC 
Ha-Hatl, dL (3) 
FX XH th Hdx + 站 Na (4) 


x 为 圆柱 形 传 递 线 ( 管 ) 内 总 的 基本 广 延 量 , E 为 圆柱 形 传 递 线 ( 管 ) 内 的 总 能 量 , ,=J,4 


为 单位 时 间 通 过 横 截面 积 4 的 基本 广 延 量 流 ，x, 为 单位 长 度 内 基本 广 延 量 的 产生 量 。 


式 (3) 表 明 ， 单 位 时 间 内 一 段 圆 柱 形 传 递 线 ( 管 ) 路 中 广 延 量 的 变化 量 ， 等 于 通过 横 截 面 
流入 该 段 圆柱 形 传递 线 ( 管 ) 路 的 广 延 量 与 广 延 量 在 这 段 圆柱 形 传递 线 〈 管 ) 路 内 流动 时 的 
产生 量 之 和 。 

式 (4) 体 现 了 不 同形 式 能 量 的 传递 共性 ， 是 研究 时 变 线 〈 管 ) 路 传递 规律 的 重要 依据 。 该 
式 中 各 项 的 物理 意义 : 左 侧 是 单位 时 间 内 圆柱 形 传递 线 〈 管 ) 路 内 总 能 量 的 变化 量 ， 右 侧 第 
一 、 第 二 项 是 单位 时 间 内 经 边界 流入 线 〈 管 ) 路 的 能 量 ， 右 侧 的 第 三 项 代表 与 强度 量 下 降 相 
联系 的 能 量 减少 量 ， 右 侧 第 四 项 代表 随 着 基本 广 延 量 的 增加 而 增加 的 能 量 。 由 该 式 可 知 ， 在 
能 量 的 传递 过 程 中 普遍 存在 着 存储 或 消散 效应 《〈 左 侧 项 )、 阻 力 及 惯性 效应 〈 右 侧 第 三 项 )。 

1) 能 量 传递 过 程 的 存储 或 消散 效应 。 在 时 变 能 量 网 络 中 ， 圆 柱 形 传递 线 〈 管 ) 内 的 总 能 
量 E 是 随时 间 变 化 的 ， 或 储存 或 消散 ， 这 些 能 量 的 存储 或 消散 效应 ， 都 是 通过 传递 介质 本 身 
来 完成 的 。 

2) 能 量 传递 的 阻力 及 惯性 效应 。 右 侧 第 三 项 ， 代 表 与 强度 量 下 降 相 联系 的 能 量 减少 量 。 
强度 量 下 降 反 映 的 是 能 量 在 传递 过 程 普 遍 受 到 的 阻力 及 惯性 效应 ， 而 与 强度 量 下 降 相 联系 的 
能 量 减 少量 ， 则 反映 了 传递 过 程 中 为 了 克服 阻力 及 惯性 而 消耗 及 转换 的 能 量 。 

2.2 时 变 能 量 网 络 的 基本 参数 
(1) 传 阻 R 

自然 界 各 种 形式 的 能 量 ， 在 传递 过 程 中 都 会 普遍 受到 一 种 阻碍 的 作用 。 这 种 阻力 效应 ， 
对 于 不 同形 式 的 能 量 有 着 不 同 的 表现 形式 ， 例 如 电力 网 络 中 的 电阻 ， 流 体 网 络 中 的 流 阻 与 及 
热力 网 络 中 的 热 阻 等 [13-15]。 能 量 在 线 ( 管 ) 路 内 传递 的 过 程 中 ， 为 克服 这 种 普遍 存在 的 阻 
力 效应 ， 将 产生 一 种 与 强度 量 下 降 相 联系 的 能 量 损耗 及 转换 。 

与 电阻 定义 相似 , 传 阻 的 定义 传递 线 ( 管 ) 路 两 端的 强度 量 差 AX 与 通过 的 广 延 量 右 之 比 ， 
即 : 


We 
H 


(5) 
(2) 传 容 C 


在 传递 过 程 中 ， 时 变 能 量 网 络 线 ( 管 ) 路 内 广 延 量 和 能 量 是 时 刻 变 化 的 。 广 延 量 及 能 量 的 存 
储 或 消散 ， 均 是 通过 传递 介质 来 实现 。 这 种 传递 过 程 中 普遍 存在 的 能 量 存 储 或 消散 效应 ， 可 


以 通过 一 个 传 容 5 的 参数 来 进行 描述 。 与 电容 定义 相似 ， 传 容 5 的 定义 为 : 


C=-H/(T) (6) 

传 容 代 表 了 广 延 量 与 能 量 在 传递 过 程 中 的 一 种 存储 或 消散 效应 ， 其 大 小 与 强度 量变 化 的 
快慢 有 关 。 
(3) 传 感 亏 


在 时 变 能 量 网 络 中 ， 由 于 惯性 的 作用 ， 任 何 发 生 高 速 瞬 态 变化 的 地 方 将 使 广 延 量 加 速 或 
减速 并 引起 强度 量变 化 。 这 种 传递 过 程 中 普遍 存在 的 惯性 效应 ， 可 以 用 一 个 传 感 工 的 参数 来 
描述 ， 其 定义 为 : 


LAX/ (7) 
传 感 是 能 量 传递 惯性 效应 的 一 种 反映 ， 其 大 小 与 广 延 量变 化 的 快慢 有 关 。 与 传 阻 反映 与 
强度 量 下 降 相 联系 的 能 量 损 失 不 同 ， 传 感 反映 的 是 与 强度 量 下 降 相 联系 的 能 量 转换 。 电 力 网 
络 中 电感 将 电场 能 转化 为 了 磁场 能 ， 流 体 网 络 流感 将 势能 〈 压 能 ) 转换 为 了 动能 。 此 外 ， 需 
要 特别 指出 ， 热 运动 能 瞬间 启动 ， 在 热能 传递 过 程 中 不 存在 惯性 效应 [16]。 
(4) 传导 G 
能 量 沿 较 长 且 密 封 性 不 好 的 线 ( 管 ) 路 传递 时 ， 有 可 能 会 出 现 泄露 ， 导 致 线 ( 管 ) 路 内 


能 量 减少 。 这 些 能 量 传递 过 程 出 现 的 泄露 效应 ， 可 以 由 一 个 传导 4 的 参数 来 描述 


3 时 变 传递 线 〈 管 ) 路 的 集中 参数 模型 
在 时 变 传递 线 〈 管 ) 路 中 ， 当 与 能 量 传递 相关 的 参数 〈 广 延 量 及 强度 量 ) 波长 4 远 远大 于 
线 〈 管 ) 路 的 特征 长 度 到 时， 这 些 参 数 的 波幅 仅仅 是 时 间 的 函数 ， 而 与 空间 的 位 置 无 关 。 这 


样 的 处 于 低频 脉动 下 的 传递 线 《〈 管 ) 路 ， 可 以 被 认为 是 具有 离散 特性 的 ， 可 以 通过 建立 其 相 
应 的 集中 参数 的 模型 来 描述 。 
3.1 普遍 化 广 延 量 传递 方程 

三 类 广 延 量 ( 基 本 广 延 量 、 能 量 和 烟 ) 的 传递 是 辩证 统一 的 ， 是 对 同一 个 传递 过 程 不 同 
层次 的 描述 。 文 献 [12] 根 据 这 三 者 的 关系 ， 结 合 传递 公理 ， 拓 展 了 能 量 传 递 公理 及 灿 传 递 公 
理 。 


为 不 失 一 般 性 ， 本 节 将 从 三 类 广 延 量 的 空间 平衡 方程 出 发 ， 通 过 分 析 及 讨论 ， 最 终 得 到 
普 裔 化 广 延 量 传 递 方 程 。 A 描述 ， 体 现 了 同一 个 传递 过 程 不 同 层次 
(基本 广 延 量 流 、 能 流 、 烟 流 ) 的 传递 共性 。 

籼 的 空间 平衡 方程 为 "71: 


(X—Xo)dx 
~ dt 


其 中 jp 为 介质 的 密度 ，x 为 单位 质量 介质 中 的 基本 广 延 量 ，X 为 强度 量 ，X, 为 强度 量 的 


VTS — Xo]+ ge (8) 
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我 态 值 ，J, 为 基本 广 延 量 的 流 密度 矢量 ，gus =J,*VX +(X 一 X0)g, 为 单位 体积 介质 中 烟 的 源 强 


下 


度 。 
对 上 式 两 端 进行 体积 积分 ， 化 简 可 得 : 
aE 
PE 00 XO -OX Ha + 0) 


[Hax + jx -XOdx, 


定义 ,= XH,，,=(X 一 X0)H, 分 别 为 单位 时 间 通 过 横 截面 积 4 的 能 流 和 烟 流 ， 联 立 式 
(3)、(4)、(9)， 即 有 : 


Ce 

7 —H,, +H,. 

dE 

i 一 已 +H,,, (10) 
dEx 

a =H,u — Hs tH 


其 中 万 ,、H,s、#, 分别 为 基本 广 延 量 流 、 能 流 、 烟 流 的 源 项 : 


HH | .gay = | 


Ho =| gadv =| Hdx + | Xdx, (11) 


Ho =| gidv =| HadX + | (X -XOdx, 
式 (10) 中 三 个 方程 ， 分 别 是 时 变 能 量 弟 线 ( 管 ) 路 的 基本 广 延 量 、 能 量 、 烟 传递 方程 。 
这 三 个 传递 方程 的 左 侧 , 分 别 代表 单位 时 间 内 基本 广 延 量 、 能 量 和 烂 的 变化 量 , 令 C,、C,,、C. 


分 别 代 表 基 本 广 延 量 流 、 能 流 及 烟 流 的 传 容 ， 根 据 传 容 的 定义 及 其 物理 意义 ， 式 (10) 可 以 转 
化 为 : 


CT Hy -HatH, 

i 
dxX 

Co dt =H, Hos | H,,, (12) 
dx 

Ca ed Pant Hoe 


上 式 中 三 类 广 延 量 传递 方程 , 在 形式 上 是 完全 一 致 的 , 是 对 同一 个 传递 过 程 不 同 层次 〈 基 
本 广 延 量 流 、 能 流 、 烟 流 ) 的 不 同 描述 。 对 一 个 具体 的 传递 过 程 ， 可 以 从 不 同 的 层面 来 进行 
研究 。 在 不 同 的 层次 上 的 具体 内 容 存 在 差别 ， 但 其 形式 始终 保持 着 一 致 。 根 据 这 种 传递 过 程 
形式 的 一 致 性 ， 普 遍 化 广 延 量 传递 方程 如 式 (13) 所 示 。 
oe (13) 
dt 
3.2 时 变 传递 线 〈 管 ) 路 的 等 效 传递 方程 
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普遍 化 广 延 量 传递 方程 , 实现 了 对 传递 过 程 的 统一 描述 , 体现 了 广 延 量 在 时 变 传递 线 管 ) 
路 内 的 变化 规律 ， 是 集中 参数 模型 建立 的 重要 依据 。 在 此 基础 上 ， 结 合 不 同形 式 能 量 的 具体 
传递 特性 ， 下 文 进一步 得 到 具有 实用 意义 的 时 变 传递 线 〈 管 ) 路 的 等 效 传递 方程 。 

(1) 电能 的 传递 过 程 

与 电能 传递 相关 的 强度 量 为 电势 ， 基 本 广 延 量 为 电荷 。 由 电荷 守恒 定律 可 知 ， 电 荷 的 源 
强度 g, =0， 结 合式 (11) 可 知 ， 与 电能 传递 相关 的 基本 广 延 量 ( 即 电荷 ) 的 源 强度 及 源 项 均 为 0。 

(2) 压 能 的 传递 过 程 

与 压 能 传递 相关 的 强度 量 为 压强 ， 基 本 广 延 量 为 体积 。 压 能 的 传递 过 程 中 ， 流 体 的 体积 
是 不 守恒 的 。 由 质量 守恒 定律 可 知 ， 质 量 的 源 强度 s =0 。 结 合 关系 式 v= wz ， 即 可 得 随 体 形 


式 的 体积 平衡 方程 : 


pO=VV (14) 


式 (1) 及 式 (14) 对 比 可 知 ,与 压 能 传递 相关 的 基本 广 延 量 (也 即 体 积 ) 的 源 强度 及 源 项 均 为 


(3) 热能 的 传递 过 程 
与 热能 传递 相关 的 强度 量 为 热力 学 温度 ， 基 本 广 延 量 为 彤 。 以 非 稳 态 热传导 为 例 ， 由 文 


献 [17] 可 知 ， 灶 的 源 强 度 ,上 0， 在 无 其 他 形式 能 量 渗入 的 条 件 下 ， 其 表达 式 可 以 简化 为 : 


=- V7 (15) 


代入 能 量 的 源 强 度 定义 式 可 得 gw =0， 即 热能 的 传递 过 程 中 ， 虽 与 热能 传递 相关 的 基本 


广 延 量 ( 焙 ) 的 源 项 不 为 0， 但 与 其 相关 的 能 量 〈 即 热量 ) 的 源 项 为 0。 

综 上 所 述 ， 对 于 电能 、 压 能 的 传递 过 程 ， 其 基本 广 延 量 的 源 项 〈 分 别 是 电荷 及 体积 ) 为 0; 
而 对 于 纯 导 热 的 过 程 ， 其 基本 广 延 量 的 源 项 ( 炉 ) 不 为 0， 但 其 能 量 的 源 项 〈 热 能 ) 为 0; 也 
即 ， 与 这 些 能 量 传递 线 〈 管 ) 路 对 应 的 广 延 量 传递 方程 ， 经 适当 的 简化 后 ， 均 可 转化 为 如 式 
(16) 所 示 的 形式 。 对 于 电能 、 压 能 的 传递 过 程 ， 广 延 量 流 为 基本 广 延 量 流 〈 即 电流 及 体积 流 
量 )， 对 于 热能 ， 广 延 量 流 为 能 流 〈 即 热流 )。 

CQX 
0 =- (16) 

能 源 有 许多 种 类 型 ， 实 际 工 程 中 能 量 的 网 络 化 传递 主要 依靠 电网 和 流体 网 络 〈 如 热 网 、 
燃气 管 网 等 )， 核 能 、 风 能 、 太 阳 能 、 水 能 等 都 必须 先 转化 为 电能 或 者 热能 再 传递 。 因 此 ， 本 
节 得 到 的 时 变 传递 线 ( 管 ) 路 的 等 效 传递 方程 ， 即 式 (16)， 是 建立 时 变 能 量 网 络 集中 参数 模 
型 的 关键 。 

此 外 ， 文 献 [12] 指 出 ， 任 何 形式 的 能 量 都 具有 质 和 量 两 重 属性 ， 能 量 的 质 是 烟 ， 烟 是 推动 
目的 过 程 进行 的 源 动力 ， 只 要 系统 的 某 一 个 强度 量 等 于 疲 态 值 ， 与 该 强度 量 相 应 的 烤 就 等 于 
0。 因 此 ， 在 时 变 传 递 线 〈 管 ) 路 的 集中 参数 模型 中 ， 强 度量 的 基准 值 〈 参 考 地 ) 应 该 为 与 该 


种 形式 能 量 相对 应 的 强度 量 寂 态 值 xX。。 
时 变 能 量 网 络 中 的 常见 能 量 〈 电 能 、 压 能 、 热 能 ) 的 传递 过 程 ， 均 可 以 通过 时 变 传递 线 


天 
AL 


〈《 管 ) 路 的 等 效 传递 方程 ， 即 式 (16) 来 进行 描述 ， 该 方程 体现 了 广 延 量 在 时 变 传递 线 〈 
路 内 的 变化 规律 。 根 据 式 (16)， 结 合 时 变 能 量 网 络 基 本 参数 的 定义 及 其 物理 意义 ， 即 可 建立 
时 变 传递 线 〈 管 ) 路 的 集中 参数 模型 ， 该 模型 如 图 2 所 示 。 


X HH H, XxX, 
oO 


O O 
Xs Xo 


图 2 时 变 传递 线 〈 管 ) 路 的 集中 参数 模型 


Fig.2 lumped parameter model of time-varying transfer line (pipe) 


4 能量 转换 设备 的 集中 参数 模型 


4.1 感应 电动 机 
感应 电动 机 的 工作 是 建立 在 电磁 感应 的 基础 之 上 ， 其 可 以 将 电能 转化 为 机 械 能 。 施 加 在 


感应 电动 机 两 端的 电压 感应 出 电机 磁场 线圈 的 速度 流 右 ,， 而 产生 的 转 矩 AXi 反馈 形成 电流 


8;， 其 转换 关系 如 式 (17) 所 示 。 这 种 转换 特性 可 以 用 回旋 器 GY(k) 的 能 量 转换 装置 来 表示 ， 其 
中 局 ,为 电动 机 转 矩 转换 系数 "1。 


(17) 


图 3 感应 电动 机 的 集中 参数 模型 


Fig.3 Lumped Pr model of induction motor 


在 建 并 感应 电动 机 的 动态 模型 时 ， 通 常 需要 考虑 电动 机 转速 差 、 轴 承 摩 擦 及 转子 摩擦 、 
转子 及 轴 的 转动 惯量 与 电机 轴 的 柔 度 等 特性 ， 以 上 特性 可 依次 由 集中 参数 元 件 R,、R,、 荆 , 


与 Cw 来 表示 "…。 由 此 ， 即 可 建立 感应 电动 机 的 集中 参数 模型 ， 该 模型 如 图 3 所 示 。 


4.2 离心 孙 
离心 录 是 流体 网 络 中 常见 的 能 量 转换 设备 ， 其 可 以 将 机 械 能 转化 为 流体 压 能 ， 其 能 量 转 
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换 特 性 能 用 回旋 器 GY(r,) 来 表示 ， 其 能 量 转换 关系 见 下 式 。 


ms (18) 


H,r, = AX, 


其 中 AX, 、 瓦 ,分 别 为 转 矩 与 角速度 ，X, 、H， 分 别 压力 与 体积 流量 ，/, 为 叶轮 回旋 器 模 数 ?0。 

在 建立 离心 泵 的 动态 模型 时 ， 通 常 需要 考虑 离心 泵 转子 摩擦 特性 、 转 子 转动 惯量 特性 、 
泄露 特性 与 阀门 出 口 压力 损失 特性 等 ,以 上 特性 可 依次 由 集中 参数 元 件 R, 、L,, 、R,, 与 R ,来 
表示 中 。 综 上 所 述 ， 离 心 泵 的 集中 参数 模型 如 图 4 所 示 。 


R L,, 
nf 1 H H 


四 qv 各 
R 
GY Ee 
AxX7 (7, )>AXn R pe 
0 oO 


图 4 离心 泵 的 集中 参数 模型 


Fig.4 Lumped parameter model of centrifugal pump 


5 时 变 能 量 网 络 方程 的 建立 


5.1 典型 时 变 能 量 网 络 

在 流体 网 络 理论 中 ， 与 基 尔 霍 夫 电 压 定律 (KVL) 和 基 尔 霍 夫 电 流 定律 (KCL) 对 应 的 分 别 是 
流体 节点 定律 及 流体 回路 定律 (时 ， 能 量 网 络 理论 中 广义 基 尔 霍 夫 定 律 的 提出 ， 使 得 我 们 能 
像 电 网 络 理论 一 样 ， 列 出 能 量 网 络 方程 中 。 为 考察 时 变 能 量 网 络 的 动态 过 程 ， 本 文 在 时 变 能 
量 网 络 模型 的 基础 上 ， 提 出 通过 构建 时 变 能 量 网 络 方程 (包括 状态 方程 和 输出 方程 》 对 多 能 
互补 系统 的 动态 特性 进行 建 模 仿真 的 分 析 方 法 ， 并 通过 算 例 分 析 对 所 提 方 法 的 有 效 性 及 实用 
性 进行 验证 。 

为 不 失 一 般 性 ， 本 文 以 某 冷 电 联 供 系统 (典型 的 时 变 能 量 网 络 ) 为 研究 对 象 ， 该 系统 如 
图 5 所 示 。 冷 电 联 供 系 统 存在 着 多 种 形式 的 能 量 传 递 ( 电 能 、 压 能 、 冷 能 等 )， 这 些 能 量 除 了 
通过 传递 线 ( 管 ) 路 进行 传递 外 ， 还 通过 多 个 能 量 转换 设备 看 合 在 一 起 (离心 隶 、 水 冷 机 组 、 
换 热 设备 等 )。 

图 5 中 实 线 代表 电力 网 络 ， 传 输电 能 ， 虚 线 代 表 流 体 网 络 ， 主 要 传输 冷 能 ， 同 时 消耗 压 
能 ;电力 网 络 与 流体 网 络 依靠 离心 泵 TWM (离心 泵 通过 感应 电动 机 带动 ) 和 制冷 机 


(P=P,+P,， 其 中 记 , 代表 水 冷 机 组 提供 的 冷 量 ，P 代表 电 空调 提供 的 冷 量 ) 耦合 在 一 起 。 


5 代表 电 源 ， 工 代表 电力 网 络 中 除了 离心 泵 及 电 空 调 外 的 其 他 负荷 ，7R 为 换 热 器 。 


在 图 5 所 示 的 时 变 能 量 网 络 中 ， 冷 能 主要 依靠 流体 网 络 进行 传递 。 在 管道 较 不 长 的 情况 
下 ， 流 体 在 管道 传递 过 程 中 冷 量 损失 很 小 ， 可 认为 流体 温度 仅 在 水 冷 机 组 及 换 热 设备 发 生冷 
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量 交换 时 才 发 生变 化 。 此 外 ， 在 实际 交流 电力 网 络 中 ， 其 电流 与 电压 通常 可 表达 为 正弦 或 者 
余弦 函数 ， 每 时 每 刻 都 在 发 生变 化 ， 但 其 作用 在 电动 机 与 电 空 调 的 电压 通常 可 以 认为 是 恒定 


的 。 为 分 析 方 便 ， 在 本 文 假定 电力 网 络 电 压 源 5, 为 常量 。 


人 TR [----<——- 


图 5 时 变 能 量 网 络 
Fig.S Time-varying energy network 
5.2 时 变 能 量 网 络 模 型 
根据 以 上 分 析 ， 结 合 网 络 相 关 知识 ， 将 图 5 简化 ， 并 将 各 个 文 路 ， 节 点 以 及 基本 回路 编 
号 ， 即 可 建立 时 变 能 量 网 络 简化 图 。 简 化 图 如 图 6 所 示 ， 其 中 子 图 (a) 为 电力 网 络 和 流体 网 络 
简化 图 ， 子 图 (b) 为 能 量 转 换 设备 (感应 电动 机 和 离心 汞 〉 简化 图 。 


Wi bs ns 


(a) 电力 网 络 和 流体 网 络 简化 图 


nn no b, mo bs mi Li 1103 
人 一 人 一 一 一 一 全 > co > oo 
= \ a a < \ ” bi 
b b 及 [多 b b 
er; 和 y < TF,) > 
{ 入 


(b) 能 量 转换 设备 简化 图 
图 6 时 变 能 量 网 络 简化 图 
Fig.6 Simplified diagram of the energy network 


在 时 变 能 量 网 络 简化 图 中 ， 选 定 n、n,、n 为 参考 节点 ， 即 可 得 到 与 之 对 应 的 关联 矩阵 4 


及 基本 回路 矩阵 B,。 
对 于 图 6， 其 广 延 量 矩 阵 为 ， 


文 路 两 端 强度 量 差 值 矩 阵 为 : 
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AX=[AX, AX, ... AX's AXio] 


全 =0 (19) 


BAX=0 


该 方程 组 中 ， 共 有 19 个 方程 ， 有 变量 38 个 ， 因 此 ， 还 差 19 个 方程 才能 解 此 方程 组 。 在 
电力 网 络 中 ， 电 压 源 5, 通常 是 已 知 的 。 在 图 5 所 示 的 时 变 能 量 网 络 中 ， 电 力 网 络 及 流体 网 络 


通过 离心 条 《感应 电动 机 ) 和 电 空 调 耦 合 在 一 起 ， 根 据 能 量 守 恒定 律 ， 离 心 条 《感应 电动 机 ) 
文 路 和 电 空 调 文 路 可 以 再 列 出 6 条 方程 。 此 外 ， 对 于 流体 网 络 、 感 应 电动 机 与 离心 泵 的 各 个 


支 路 ， 可 再 建立 12 条 支 路 特性 方程 。 流 体 网 络 支 路 对 应 的 传 阻 R,( 层 流 状态 ) 与 R,，( 闲 


流 状态 )、 传 容 C 及 传 感 工 可 由 式 (20) - (22) 求 得 。 联 立 上 述 方程 ， 即 可 建立 时 变 能 量 网 络 
模型 。 


_ 1284L 
lam 4 
8 人 
p 
RR, 局 元 2 站 5 H,, 


AL 
Ca (21) 


pL 
RY ed 22 
2 2) 
上 3 式 中 ， 为 动力 粘度 ，L 与 D 为 管道 长 度 及 直径 ;pp 为 流体 密度 ，H， 为 流 过 管道 的 


体积 流量 ，4 为 率 流 状态 下 的 管 路 沿 程 阻力 系数 ， 其 能 采用 Karman-Nikuradse 光滑 管 公式 进 
行 计算 ”A 为 管道 横 截面 积 ，a 为 流体 的 扰动 波 传播 速度 。 
5.3 时 变 能 量 网 络 方程 

系统 在 输入 或 扰动 作用 下 的 运动 状态 可 用 一 组 状态 变量 来 描述 ， 状 态 变量 是 借以 表征 系 
统 内 部 状态 随时 间 变 化 的 物理 变量 ， 而 描述 状态 变量 随时 间 变 化 的 数学 表达 式 即 状态 空间 方 
程 。 为 考察 时 变 能 量 网 络 的 动态 过 程 ， 本 文 在 时 变 能 量 网 络 模型 的 基础 上 ， 进 一 步 构建 时 变 
能 量 网 络 方程 (包括 状态 方程 与 输出 方程 ) 并 利用 数值 仿真 软件 对 冷 电 联 产 系统 进行 动态 仿 
真 及 分 析 。 

(1) 控制 变量 x 与 状态 变量 x* 的 选取 

对 于 如 图 6 所 示 的 时 变 能 量 网 络 ， 施 加 在 电动 机 与 电 空 调 两 端的 电压 源 5, 是 外 加 的 激励 


源 ， 因 此 控制 量 w = AX,。 


一 般 来 说 ， 状 态 变量 选取 的 个 数 通常 与 系统 储 能 元 件 的 个 数 一 致 。 但 如 果 系统 某 些 储 能 
元 件 依赖 于 其 他 储 能 元 件 时 ， 则 状态 变量 数 要 少 于 储 能 元 件数 。 


对 如 图 6 所 示 的 时 变 能 量 网 络 ， 可 分 别 选 取 传 感 文 路 的 广 延 量 值 8;、H6、H;、H'is、 H's 
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与 传 容 支 路 的 强度 量 差 和 AX,。、AX。、AX,,。、 AX,, 作为 状态 变量 。 


(2) 根据 4.2 节 建 立 的 时 变 能 量 网 络 方程 ， 进 一 


方程 及 输出 方程 。 


特别 地 ， 对 于 线性 系统 ， 方 程 (23) 可 进 一 


= f (x,u,t) 
y= 8(X,u,1) 


X= Ax+Bu 
y=Cx+Du 


呈 些 元 素 是 时 间 的 函数 


步 推导 出 如 式 (23) 的 时 变 能 量 网 络 状 态 


步 写 成 方程 (24) 所 示 的 标准 形式 。 


(24) 


， 式 (24) 的 系数 矩阵 的 各 元 素 丝 为 常数 ， 吞 系统 为 时 变 系统 ， 则 系数 矩阵 中 


6.1 时 变 能 量 网 络 的 基本 参数 
仍 以 上 述 冷 电 联 供 系 统 为 例 ， 该 系统 总 的 冷 负 和 荷 为 B.=500kW ， 该 冷 负 荷 由 水 冷 机 组 和 中 
央 电 空调 系统 共同 提供 。 该 系统 能 量 转换 设备 (异步 感应 电动 机 、 离 心 泵 及 电 空调 ) 的 基本 
参数 如 表 1 所 示 。 
表 1 能 量 转 换 设备 的 基本 参数 
Tab.1 Parameters 
元 件 运行 参数 单位 描述 
R=53 N:m/(rad l/s) 滑 差 系 数 
KRorl=3.16x101 Nm/ (rad /sy 旋转 摩擦 系数 
步 由 到 
异步 电动 机 KR 32x10-4 N:m/(rad/s) 东 摩擦 系数 
Lnx =1. 65x10 kg :m2 转子 惯量 
Cr =103 rad / (Nm) 轴 和 柔 度 
Ror =8x102 N:m/ (rad /'s) 转子 摩擦 系数 
Si Lyx =1.15x103 kg-m? 转子 惯量 
Rye =5x10!! Pa/ (m/s) 泄露 系数 
Rp, =15 阀门 阻力 系数 
空调 COP=4 能 效 系数 


6.2 仿真 分 析 


根据 上 一 章 的 研究 内 容 ， 即 可 建立 时 变 能 量 网 络 的 动态 模型 。 


本 文 在 上 述 模型 的 基础 上 ， 
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对 三 种 不 同 扰动 发 生 后 的 能 量 网 络 响应 进行 时 域 仿 真 ， 该 三 种 扰动 分 别 为 : (1) 电压 源 8. 突 
变 ， 其 电压 有 效 值 由 220V 下 降 至 218V; (2) 增 大 离心 泵 的 出 口 阀 门 开 度 ， 阀 门 阻力 系数 突 
变 为 R, =10; 〈3) 减 小 离心 泵 的 出 口 阀门 开 度 ， 阀 门 阻力 系数 突变 为 R=20。 

本 文 建立 的 时 变 能 量 网 络 方程 〈 状 态 方程 和 输出 方程 )， 由 于 考虑 了 流体 网 络 与 能 量 转换 
设备 的 非 线性 动态 特性 ， 实 质 为 一 个 非 线性 微分 方程 组 。 对 于 形 如 式 (23)、 式 (24) 的 微分 方程 
组 ， 文 献 [22] 与 [23] 提 供 了 许多 行 之 有 效 的 解法 。 本 文 利用 变 步 长 的 龙 格 库 塔 法 进行 求解 ， 并 
通过 MATLAB 进行 仿真 ， 仿 真 时 长 为 2s8。 通 过 把 异步 感应 电动 机 、 离 心 汞 、 电 空调 的 基本 
参数 以 及 初始 条 件 代 入 后 ， 经 过 和夫 代 即 可 求 得 时 变 能 量 网 络 的 动态 响应 。 

三 种 不 同 扰动 发 生 后 电力 网 络 、 流 体 网 络 与 感应 电动 机 的 动态 响应 分 别 如 图 7、 图 8 与 图 
9 所 示 。 图 7 为 电力 网 络 离心 泵 支 路 电流 已 与 电 空 调 支 路 电流 已 的 动态 响应 ， 图 8 为 流体 网 


络 离心 稍 支 路 体积 流量 尽 与 出 口 压力 AXss 的 动态 响应 ， 图 9 为 感应 电动 机 的 转速 已 , 与 转 矩 


AXi 的 动态 啊 应 。 


Tr 


图 7 电力 网 络 动 态 响 应 


Fig.7 Dynamic response of electric power network 


图 8 流体 网 络 动态 响应 


Fig.8 Dynamic response of fluid network 


INNAcreoerer er rrr 


7 一 一 一 一 一 一 一 一 一 一 一 一 一 一 


Da 


9 感应 电动 机 动态 响应 


Fig.9 Dynamic response of induction motor 


扰动 1 发 生 时 ， 在 :=0 时 刻 电压 源 5, 从 220V 突变 为 218V ， 电 压 水 平 下 降 ， 导 致 离心 泵 支 


路 电流 豆 、 离 心 泵 出 口 压力 Axis 、 电 动机 的 转 矩 AX 与 转速 已 ,及 均 发 生 大 幅度 的 脉动 。 扰 


动 1 发 生 后 ， 电 压 突 变 首 先 冲 击 感应 电动 机 ， 其 转速 已, 由 155.8rad/s 突变 至 153.5rad/s ， 其 输出 
转 抢 也 由 扰动 前 的 10.IN. 关 下降 至 8.6N.m 。 电 动机 转速 与 转 符 的 大 幅度 下 降 导 致 离心 泵 出 口 
压力 AXs 与 体积 流量 五 , 减 小 ， 扰 动 1 发 生 后 电动 机 与 离心 条 输出 功率 均 变 小 。 供 水 量 的 减少 


导致 供 冷 量 大 幅 下 降 ， 在 电压 下 降 的 情况 下 ， 为 满足 冷 负 荷 需求 ， 电 空调 文 路 电流 太 , 出 现 较 


大 幅度 增加 。 
扰动 2 发 生 时 ， 离 心 泵 出 口 阀门 开 度 突然 增 大 ， 阀 门 对 离心 人 汞 产生 的 阻力 减 小 ， 离 心 泵 


体积 流量 五 , 增 大 。 根 据 离 心 录 的 特性 曲线 可 知 ， 通 过 离心 泵 的 流量 增 大 ， 离 心 泵 输出 的 有 效 
扬程 将 减 小 ， 该 特性 与 图 8 流体 网 络 动态 响应 相 一 致 ”。 流 体 网 络 运行 工 况 的 变化 ， 将 反 过 
来 影响 电力 网 络 和 感应 电动 机 参数 。 扰 动 2 发 生 后 ， 离 心 泵 支 路 电流 五 ,和 感应 电动 机 有 效 转 


矩 AXi, 发生 持 续 下 降 ， 分 别 由 初始 时 刻 的 8.04 和 10.1N.m 下 降 至 7.74 和 9.7N.m。 此 外 ， 由 于 


离心 泵 供水 量 不 断 瑞 升 ， 电 制冷 需求 减少 ， 电 空调 支 路 电流 已 发 生 持续 性 下 降 。 
扰动 3 发 生 时 离心 泵 出 口 阀门 开 度 突然 减 小 ， 与 扰动 2 相反 ， 阀 门 对 离心 泵 产生 的 阻力 
增加 ， 离 心 泵 体积 流量 总 减 小 出 口 压力 AXis 增 大 。 扰 动 2 发 生 后 ， 离 心 泵 文 路 电流 已 和 感 


应 电动 机 有 效 转 矩 AXi, 持续 上 升 ， 分 别 由 初始 时 刻 的 8.04 和 10.1N.m 增 大 至 8.34 和 10.5N.m。 


由 图 9 可 知 ， 当 扰动 2>、3 发 生 时 ， 由 于 电压 源 5, 恒 定 在 220V ， 流 体 网 络 对 电动 机 转速 石 ;的 
影响 很 小 ， 可 近似 认为 转速 恒定 。 此 外 ， 由 于 离心 泵 供水 量变 少 ， 电 制冷 需求 增 大 ， 电 空调 
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文 路 电流 已 发生 持续 性 上 升 。 


由 此 证 明了 本 文 所 提 理 论 的 有 效 性 及 实用 性 。 本 文 提 出 的 时 变 能 量 网 络 动态 建 模 及 仿真 方法 ， 


由 上 述 的 分 析 可 知 ， 基 于 本 文 所 提 方 法 的 仿真 结果 与 冷 电 联 产 系统 的 运行 特性 是 一 致 的 ， 


对 多 能 互补 系统 的 优化 设计 及 性 能 分 析 具 有 重要 的 实际 意义 。 
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结论 
在 已 有 能 量 网 络 理论 基础 上 ， 从 能 量 本 质 的 角度 出 发 ， 对 时 变 传递 线 〈 管 ) 路 和 能 量 耦 


合 设备 的 能 量 传递 机 理 进 行 深 入 探讨 ， 建 立 了 时 变 能 量 网 络 模型 。 为 考察 时 变 能 量 网 络 的 动 
态 过 程 ， 本 文 在 时 变 能 量 网 络 模 型 的 基础 上 ， 进 一 步 构建 时 变 能 量 网 络 方程 (包括 状态 方程 
和 输出 方程 )。 基 于 本 文 所 提 方 法 对 茶 冷 电 联 产 系统 在 不 同 扰动 发 生 后 的 动态 啊 应 进行 数值 仿 


真 ， 仿 真 结 果 与 冷 电 联 产 系统 的 动态 特性 一 致 ， 由 此 证 明了 本 文 所 提 分 析 方 法 的 有 效 性 及 实 
用 性 。 本 文 的 研究 内 容 为 时 变 能 量 网 络 的 建 模 、 分 析 、 优 化 及 规划 莫 定 了 基础 。 
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